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ABSTRACT 

The recent detection of the anomalous X-ray pulsar (AXP) 4U 0142+61 in the 
mid- infrared with the Spitzer observatory by Z.Wang and coworkers constitutes 
the first instance of a disk around an AXP. We show, by analyzing earlier optical 
and near-IR data together with the recent data, that the overall broadband data 
set can be reproduced by a single model of an irradiated and viscously heated 
disk. 

Subject headings: pulsars: individual (AXPs) — stars: neutron — X-rays: bursts 
— accretion, accretion disks 



INTRODUCTION 



The zoo of young neutron stars contains a number of categories recognized by the 
distinct properties of such stars discovered or identified in the last decade. The existence 
of these different types of neutron stars strongly suggests that not all neutron stars are 
born with the same initial conditions, nor do they follow the same evolutionary path as 
the familiar radio pulsars, of which the typical young example is usually taken to be the 
Crab pulsar. Anoma l ous X -ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) (see 
Woods &: Thompson 2005!. for a review of AXPs and SGRs) are widely accepted to be 



magnet ars ( 



Duncan fc Thompson! 1 19921 ). The "radio-quiet" neutron stars (also known as 



central compact object s [CCOs]; (P avlov et al. 2004) and dim isolated (thermal) neutron 
stars (DTNs or DINs; iHaberll |2005| ) are probably related classes, although the latter are 
relatively older. M easured rotat ion periods of AXPs, SGRs and DINs cluster in the narrow 
range of 3 — 12 s. lAlparl (1200 ll ) proposed that the presence or absence, and properties, of 
bound matter with angular momentum that is, a fallback disk, may be among the initial 
parameters of newborn neutron stars, in addition to magnetic dipole moment and initial 
rotation rate. According to this scenario, the differences between isolated pulsars, AXPs, 
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SGRs, and DINs, as well as the CCOs in certain supernova remnants, are due to different 
initial conditions, including the presence or absence and properties of fallback disks. 

The suggestion that the X-ray luminosity of AXPs is due to mass accretion from a 
fallback disk (Chatterjee et al. 2000; Alpar 2001) has triggered searches for evidence of 
such disks. Searches for fallback disks around AXPs have been conducted primarily in the 
optical and near IR bands. A realistic model for a putative fallback disk should take into 
account the effects of irradiation by the neutron star's X-ray luminosity. Irradiation is the 
dominant source of disk luminosity for the outer disk at IR and longer wavelengths. The 
radial temperature profile of a fallback disk has been co mputed by Pe r na et al. (2000) and 



HuUeman et al. (2000) using the prescription given by IVrtilek et al.l (Il990l ). To estimate 
the temperatures, both HuUeman et al. and Perna et al. assumed the same particular 
irradiation efficiency and found that their estimated optical flux values lie well above the 
values indicated by observations of the AXPs 4U 0142+61 and IE 2259+586. While Perna 
et al. suggested that the difference might be due to a probable advection-dominated flow at 
the inner disk, HuUeman et al. drew the conclusion that this model can only fit the data for 
an unrealistically small outer disk radius and that fallback-disk models were ruled out (see 
§ 4 for further discussion). Recent studies show that irradiated-disk model fits to existing 
detections and upper limits from AXPs allow th e presence of accreting fallback disks with 



reasonable irradiation parameters for all AXPs (lErtan fc Cali§kanl |2006| ) when the fits are 
not restricted to a particular irradiation efficiency. 

Recent results by Wang et al. (2006, hereafter WCK06) from Spitzer data on the AXP 
4U 0142+61 complement earlier data in the near-infrared and optical (HuUeman et al. 2000, 
2004, hereafter HvKK OO and HvKK04) and the detection of pulsed emission in the optical 



( iKern fc Martini l2002l ). thus providing the first instance of a multiband data set that can be 
tested against disk models. WCK06 interpret their data with a dust disk radiating in the 
infrared, while the optical data are ascribed a magnetospheric origin. They take the origin 
of the X-ray luminosity to be dissipation of magnetic energy in the neutron star's crust, in 
accordance with the standard magnetar model. 

Here we show that the entire unpulsed optical, near-IR, and far-IR spectrum can be 
explained as due to a gaseous disk whose luminosity is provided by viscous energy dissipation 
due to mass transfer through the disk, as well as irradiation by the X-ray luminosity Lx from 
the neutron star. Taking Lx to be an accretion luminosity, as we do, links the irradiation 
contribution to the disk luminosity with the mass accretion rate Mace and, thereby, the mass 
inflow rate M through the disk and the viscous dissipation contribution of viscous dissipation 
to the disk luminosity. Luminosity in the optical band comes from the inner regions of the 
disk and is supplied, to a significant extent, although not fuUy, by the local viscous energy 
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dissipation in the inner disk. The innermost disk, which is dominated by dissipation, radiates 
mostly in the UV band but also contributes to the total optical emission. Extension of the 
same irradiation model to outer regions of the disk explains the IR radiation. The disk 
luminosity is dominated by the contribution of irradiation at outer radii, where the thermal 
radiation peaks in the longer wavelength bands (see Table 1). 



Pulsed emission in the optical (IKern fc Martirul2002l ) is due to mag netospheric emission 
power ed by a disk-shell dynamo in our model. It has been shown by I Cheng fc Ruderman 
( 1l99l[ ) that magnetospheric pulsar emission in optical and higher energy bands can be sus- 
tained by an outer gap in the presence of a disk protruding int o the magnetosphere. T his disk 
model is different from the first disk- magnetosphere models (IMichel &: Desslerlll98ll ). More 
recently, lErtan &: Chengl (120041 ) showed that the pulsed optical emission from 4U 0142+61 
can be explained in terms of a magnetospheric model with a disk and a neutron star surface 
dipole magnetic field of -Bdipoic ~ 10^^-10^'^ G. The magnetospheric emission is nearly 100% 
pulsed in this magnet osphere-disk model. The discovery of strongly pulsed optical emission 
( iKern fc Martinll2002l ) has generally been perceived as excluding disk models, on the basis of 
the notion that magnetospheric emission cai inot survive with a di s k pro truding into the light 
cyl inder. The disk- r nagne tosphere model of ICheng fc Ruderman! (jl99ll ) and the application 
by lErtan fc Chengl (120041 ) are important for fallback-disk models, since they demonstrate 
that a disk can actually be part of a magnetospheric mechanism that generates strongly 
pulsed emission. 

The clustering o f AXP and SGR rotational periods between 6 and 12 s was investigated 
by iPsaltis &: Millen (120021 ). who found that statistical clustering with a spin-down law of 
braking index n ~ 2-4 implies that the final period of these systems must be about 12 s. This 
final period does not correspond to reaching the "death line" in the P-P plane which refiects 
the cut-off voltage for magnetospheric emission. Within isolated dipole-magnetar models, the 
alternative possibility for reaching a final clustering period, asympt otically, is by magneti c 
field decay on timescales commensurate with the AXP and SGR ages (IPsaltis fc Millerll2002l ). 
An investigation of various magnetar field decay models (Colpi et al. 2000) has shown that 
having the X-ray luminosity data be associated with magnetic field decay is consistent with 
only one model, involving crustal Hall cascade, and for decay times less than about 10^ yr. 
The pr esence of a fallback disk prov ides a natural explanation for period clustering of the 
AXPs ( IChatterjee fc Hernguistl 120001) . as well as dim isolated (thermal) neutron stars (see, 
e.g., Haberl 2005) and SGRs JAlpailboOll : lEksi fc Alpai1l2003h . Fallback-disk models do not 
address the origin of bursts, while magnetar models explain the bursts by employing the 
energy stored in very strong magnetic fields. The basic workings of the magnetar models 
require only local magnetic stresses and energies in the neutron star crust, and they can thus 
be sustained with higher multipole components of the field, while the period clustering and 
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explanation of torques with a fallback disk require the dipole conaponent of the magnetic field 



to be on the order of 10^^-10^^ G. Hybrid models (jEk§i &: Alpar 



20031 : lErtan fc Alpai]l2003f ) 



with -Bmultipole 

fallback disk. 



10^^-10^5 G and B, 



dipolc 



10^^-10^^ G are a possibility that incorporates a 



In this paper, we evaluate the recent observations by WCK06 in conjunction with older 
results (§ 2). We show that the data agree quite well with a gaseous disk. In § 3, we estimate 
the characteristic lifetime of such a disk that has evolved through a propeller phase and find 
that the disk's estimated age, depending on the torque model, is comparable to the P^/2P^ 
characteristic age estimate taken as if dipole spin-down were effective. The discussion and 
conclusions are given in § 4. 



2. INFRARED AND OPTICAL RADIATION FROM THE DISK AROUND 

4U 0142+61 

Recent Spitzer observations of the anomalous X-ray pulsar 4U 0142+61 have detected 
the source in the 4.5 and 8 /xm infrared bands (WCK06). We will show that, taken to- 
gether with the earlier observations in the optical (HvKKOO, HvKK04; Morii et al. 2005), 
these detections are consistent with the presence of a gaseous disk around the neutron star. 
The optical detection must have a contribution from the magnetosphere, as implied by the 
strongly pulsed part of the optical signa l , at 2 7% amplitude in the R band, observed by 



Kern fc Martini (120021 ). lErtan fc Chengi (12004J ) showed that the strongly pulsed emission 
from the magnetosphere can be accounted for either with a dipole-magnetar-based magneto- 
sphere with surface dipole field strength -Bdipoie ~ 10^^-10^^ G or with a magnetosphere with 
a disk in it, with -Bdipoie ~ 10^^-10^^ G. The disk-magnetosphere model can provide nearly 
100% pulsed magnetospheric emission in the optical. The presence of pulsed magnet ospheric 
emission does not rule out the presence of a disk protruding into the magnetosphere. On the 
contrary, such a disk can be part of t he magnetospheric circuit sustai ning pulsed emis s ion, a s 
shown bv lCheng fc RudermaJ Jl99lh . and applied to 4U 0142+61 bv lErtan fc Chend kooi ). 
A pulsed component is likely to be present in the B, V, and IR bands. Uncertainties due to 
an unresolved pulsed fraction that is constant across the bands cannot be distinguished from 
uncertainties in other overall factors such as distance and disk inclination angle. Possible 
variations in the relative pulsed fractions in different bands should be taken into account 
in evaluating the model fits. Further, an important uncertainty arises from the interaction 
between inner disk and magnetosphere, which might distort the inner disk's geometry and 
thus possibly lead to time- dependent self-shielding at some inner disk regions. Observations 
of various AXPs in some of the optical/IR bands at different epochs show considerable same- 
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source variations in relative amplitudes (see, e.g., Morii et al. 2005). These variations are 
not included in our or other scientists' model fits, and indeed, it is not feasible to model such 
effects in a plausible way. Keeping this in mind, we try to obtain the best fit to optical and 
IR data points. 

Figure 1 shows the best model fit to all the data. Table 1 details the features of this 
fit. For each radiation band, the effective blackbody tcmpcratiirc representing that band 
and the radius R in the disk corresponding to that effective temperature are listed. The 
last column gives the ratio of dissipation- to irradiation-powered disk flux at that radius. 
The total emission in each band is predominantly supplied by parts of the disk near the 
quoted radius, but of course, there are lesser contributions from all parts of the disk. In 
the best fit displayed in Figure 1 and Table 1, the representative radius, 1 x 10^0 cm, that 
dominates the V-band emission draws about 30% of its luminosity from viscous dissipation. 
The actual inner disk radius in this fit. Tin = 10^ cm, is almost an order of magnitude smaller 
than the V-band representative radius. In the inner disk, the ratio of viscous dissipation to 
irradiation-supported flux is therefore higher. The radiation from the innermost disk is a 
major contributor to the UV radiation bands. The innermost disk, from Tin to the V-band 
representative radius, also makes important contributions to the radiation in the optical 
bands. Therefore, the fits are quite sensitive to the value of the inner disk radius. The same 
disk model accommodates the earher optical and IR data (HvKKOO, HvKK04; Morii et al. 
2005) and the Spitzer detections by WCK06. The value of Ay indicated by the best model 
fit shown in Figure 1 is 3.5. This is within the reasonable dereddening range 2.6 < Ay < 5.1 
for 4U 0142-F61 (HvKKOO, HvKK04) and coincides, fortuitously, with a recently derived 
estimate of Ay = 3.5 ± 0.4 (Durant & van Kerkwijk 2006). For this model, all the model 
flux values remain within ~ 30% of the data points. The error bars on the B-band data 
point give the 3 a upper limit (HvKK04). Considering the uncertainties discussed above, 
the irradiated gaseous disk model is in reasonable agreement with optical and IR data. 

In the inner disk, viscous dissipation contributes significantly to the disk luminosity, 
while in the outer and cooler parts of the disk, the luminosity is determined by X-ray 
irradiation from the neutron star. Notably, the total disk luminosity of a gaseous disk, 
including both irradiation and the intrinsic dissipation, does not fall off in the optical range 
in the way expected for a dust disk as employed by WCK06 in their Figure 3. Indeed, a 
gas disk model can explain both the optical and the IR data, as shown in Figure 1. The 
earlier conclusion by HvKKOO that an accretion disk is ruled out is no longer valid given 
the Spitzer observations. These authors ruled out an accretion disk because their optical 
magnitude, dereddened with Ay = 5.4, does not fall on the disk curve corresponding to 
Lx as an accretion luminosity, with certain assumptions, a choice of M, and a particular 
irradiation reprocessing efficiency. Fitting the optical points to the Rayleigh- Jeans end of a 
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disk blackbody therefore required, on the basis of the (only optical and near-IR) data then 
existing, that any disk in 4U 0142+61 be truncated at a certain outer radius rout- Now, with 
the new Spitzer observations of WCK06, the disk is actually detected in the far-infrared. 
We can fit all the data in the optical and IR bands B, V, R, I, J, H, K^, 4.5 /xm, and 8 fim 
(Fig. 1). A gaseous disk emitting in this energy range extends from an inner disk radius 
Tin — 1 X 10^ cm to an outer radius of about 10^^ cm (see Table 1). The two mid-IR Spitzer 
data points put a lower limit on the outer disk radius of around 10^^ cm, but they do not 
constrain the extension of the outer disk beyond 10^^ cm. In our model, the outer disk radius 
is Tout = 2 X 10^^ cm. The quality of the fits does not change for even larger values rout, since 
the contribution to the Spitzer bands at 4.5 and 8 fim from larger radii is not significant. 

The model fits in the optical band are sensitive to the value of the disk's inner radius. 
This is taken to be the Alfven radius, 

rin = rA = hM-'/\GM,)-'/'fif, (1) 

where is the mass and yU* is the dipole magnetic moment of the neutron star. The factor 
/i ~ 0.5-1 describes the uncertainty in the precise location of the inner radius. The mass 
inflow rate M in the disk is larger than or equal to the mass accretion rate derived from the 
X-ray luminosity, 

M_ ^ ^ ^ f,M < M, (2) 

where i?* is the radius of the neutron star. As the star is spinning down while accreting, it 
must be in the propeller re gime but clo se to rotational equilibrium with the disk, according 



to the fallback-disk model (lAlparll200ll ). Some of the mass inflow through the disk may not 
be accreting under these circumstances. From the observed X-ray luminosity of Lx — 10^^ 
ergs s~^, taking = 10^ cm and = 1.4 we obtain Mace — 5 x 10^^ g s~^, which is 
related to the mass inflow rate M, depending on the choice of f2- We first report model fits 
with M = Mace if 2 = 1) and discuss modifications of the results for (/i, < 1 in § 4. For 
a given M, the inner radius depends on the strength of the dipole magnetic moment of the 
neutron star. 

We calculate the radial effective temperature profile Teff(r) considering both the disk 
surface flux due to the intrinsic viscous dissipation rate D of the disk, 

D^±^ (3) 

87r r-^ 

(see, e.g., Frank et al. 2002) and the irradiation flux. 



M 
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( IShakura &: SunyaevI Il973l ). where the irradiation parameter C includes the effects of disk 
geometry, albedo of the disk surface, and the conversion efficiency of the mass accretion into 
X-rays. In our calculations, we take C to be constant along the disk and leave its value as 
a free parameter. For comparison with the data, we i ntegrate the disk blackb ody emission 



in each of the optical and IR observational bands (see lErtan fc Caliskarul2006l . for details of 
calculations). The model energy flux values presented in Figures 1 and 2 were obtained with 
C = 1 X 10~^. While the irradiation geometry of AXP disks might be somewhat different 
from those of low-mass X-ray binary (LMXB) disks (Ertan et al. 2006; Ertan & Qali§kan 
2006), the values of C that we find in our fits lie in the range expected for LMXBs (~ 10"*^ 
to 10^'^), in particular, based on disk stability analyses of soft X-ray transients (deJong et 
al. 1996; Dubus et al. 1999; Ertan & Alpar 2002). With C's in this range, the optical/IR 
observations for the other three AXPs with available data can also be fitted with disk models 



with irradiation and viscous dissipation (lErtan &: Cali§karul2006l ). 



In Figures 1 and 2, we present two different model fits. Figure 1 shows the best fit, 
obtained with Ay = 3.5, with the dipole field at the neutron star's surface set to be 10^^ G 
on the equator (2 x 10^^ G at the poles). After these fits were made, we learned of work 
on X-ray and optical extinction and reddening (Durant & van Kerkwijk 2006) reporting 
Ay = 3.5 ± 0.4 for 4U 0142-1-61. The inner disk radius r^^ for the model in Figure 1 is 10^ 
cm. The model fits are comparable for values of Tin up to a few times 10^ cm, with smaller 
Tin corresponding to larger Ay- The corotation radius for 4U 0142-1-61 is — 7 x 10*^ cm, 
which is near the inner disk radius we obtain from the best model fit. This implies that 
the system is indeed near rotational equilibrium. The system is in the propeller regime, as 
indicated by the spindown of the AXP. Nevertheless, most or all of the inflowing disk matter 
is being accreted onto the neutron star rather than being effectively propelled out of the 
system, on account of the closeness to equilibrium. This is consistent with our model flts 
obtained with the assumption that /2 = 1. 

The largest value of the dipole magnetic fleld allowed by the model fits corresponds to 
the lowest value of Ay- For the direction of 4U 0142-1-61, reasonable values of the dereddening 
in the direction of the AXP he in the range 2.6 < Ay < 5.1 (HvKKOO, HvKK04). Figure 2 
shows the model fit for the largest value of the magnetic dipole field. This was obtained 
by setting Ay = 2.6, the minimum value in the range considered, and allowing a 60% 
discrepancy in the V band. Note that the discrepancy in V is ~ 30% for the best fit, with 
Ay = 3.5. The inner radius of the disk is = 8 x lO^cm. The corresponding dipole 
magnetic field value is found to be 2 x 10^^ G on the equator (4 x 10^^ G at the poles). 
The disk model fits preclude higher values of the dipole component of the magnetic field. A 
stronger dipole field, > 10^'^ G, would cut the disk at much larger radii, above 10^° cm. 
This is not consistent with the optical (R, V, B) data. 
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In comparing dereddened fluxes with models for various values of -Bdipoie and Ay, we 
find that the highest discrepancies always occur in the V band, where the dereddened fluxes 
are always somewhat larger than the model predictions. In the B band, the observations 
have yielded only an upper limit. Predictions in the B band from models that fit the V, R 
and IR bands are consistent with this upper limit. 



3. LIFETIME OF THE DISK AROUND 4U 0142+61 

Could a fallback disk remain active, exerting a torque on the neutron star, at the age of 
4U 0142+61? In this section, we explore the evolution of the disk under various models for 
the torque between the disk and the neutron star. The characteristic lifetime of a fallback 
disk around a propeller can be estimated as = M^/M, where M is the mass inflow rate in 
the disk. For a steady configuration, the disk's mass-loss rate is Md = —M. The mass-loss 
rate will include the accretion rate Mace onto the neutron star, as well as Mout, the mass 
lost from the disk— neutron star system. In an effective propeller situation, Mout may make 
up most of Md- As the outer radius of the disk is not constrained by the data, the disk mass 
cannot be calculated from the fits. Instead, we shall assume plausible values, consistent with 
the observed disk to within an order of magnitude. We then use torque models to estimate 
the mass inflow rate M through the disk from the observed spin-down rate Cl^,. 

The long-term spin history of a rapidly rotating magnetized neutron star depends on 
the types of torques acting on it. The net torque on a magnetized neutron star interacting 
with an accretion disk can be expressed in general as 

A, = hQ, = jM{GM,r,,f'^, (5) 

where j is the nondimensional torque as a function of the fastness parameter cu* = r2*/i7K('"in) , 
n* being the angular spin frequency of the neutron star. Tin the inner disk radius, VLy. = 
{GM^/r^y^"^ the Keplerian angular velocity, and /* the moment of inertia of the neutron 
star. Using the observed value of in equation ([5]), we estimate the mass inflow rate as 

where Mn is the mass inflow rate in units of 10^'' g s~^, 1*45 is the neutron star's moment 
of inertia in units of 10^^ g cm^, and = 27r/fi* is the spin period of the neutron star 



(lErkut fc Alpai]l2004J l 



In earlier epochs of its evolution, an AXP with a fallback disk is likely to be in the 
propeller phase, in which the neutron star is rotating rapidly with respect to the rotation 
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rates in the inner boundary of the disk. Most or all of the mass supplied by the disk is lost 
from the system rather than being accreted by the propeller neutron star. This goes on until 
the neutron star has spun down to rotation rates close to corotation with the disk. When 
the star rotates only a little faster than the inner disk, a significant fraction of the mass 
lost by the disk may be accreted onto the neutron star while, at the same time, the star is 
spinning down under torques applied by the disk. 



The propeller-type spin-down torqu e first consid ered by llllarionov fc SunyaevI (Il975l ) 



scales as |A^i| ~ (/i^/r^)f2K(r'A)/f^* (see lGhoshlll995h . Depending on the spin rate of the 



neutron star as compared with the sound speed at the magnetospheric boundary, the scaling 
of the braking torques for subsonic and super sonic propellers are given by [iVol ~ u lfll/GM^, 



and |A^3| ~ iJ,in^nKirA)/GM^, respectively fiDavies fc Pringlelll98ll : lGhoshlll995h . In these 



models, the magnetospheric radius is the Alfven radius (see eq. [1]). The spin-down 
torque for a subsonic propeller, N2, is independent of M and not used in the present 
work. The torques Ni and can be expressed as in equation ([5]) with j = —u^ where 
q = —1 and q = 1 for A^^i and N^, respectively. The braking torque considered by Wang 
& Robertson (A^4; 1985) has the same form as N^, but with a particular estimate for the 
inner radius of the disk in the propeller stage, that is, N4 = —uj^M{GM^,r^^^/'^ ^ where 
= 2~^/^^VU^^^ M~^/^{2GM^Y/^ . This scaling of the inner disk radius is simply that 
of the Alfven radius when r^^ = rp^ = Vco = {GM^/illY^^, the corotation radius. All the 
braking torques yield the same M when tu* = 1, in other words, when rjn = Vc 



CO- 



We now estimate the lifetime of a fallback disk around 4U 0142+61 using these different 
types of propeller torque for disks of mass Md = 10~^ Mq and = 10~^ Mq. For each 
torque model, solving equation with the implicit M- and 5*-dependence of cj* and j{uj^) 
gives the current value of the mass inflow rate M. For = 10~^ Mq, the inferred values 
for the lifetime of the disk around 4U 0142+61 corresponding to the propeller torques A^i, 
N3, and A^4 can be written for = 10^ cm and = 1.4 as n ~ 1.0 x lO^B^^^^yi, 
Ts ~ 4.8 X lO^B^dyT, and r4 ~ 1.8 x lO^BKly^, respectively, where i?*i2 is the surface dipole 
magnetic field strength on the poles of the neutron star in units of 10^^ G. For Md = 10~^ Mq, 
all age estimates are greater by a factor of 10. In the initial life of the neutron star, after 
a brief initial accretion phase the propeller regime will start when > r^o and will prevail 
until the inner radius reaches the corotation radius again. As the inner radius of the disk r^^ 
approaches Tco, the propeller regime will allow accretion. For all torque models, the condition 
'"in = yields the same mass inflow rate, 

M = ^]p(GM,)'^/V* (7) 
which yields for = 10^^ Mq the critical age estimate 

r, 4.67 X 10^ yr B;,l (8) 
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As shown in Figure 3, Tc is the minimum estimated hfetime for a fallback disk if the neutron 
star is spinning down under the action of propeller torques without accreting any matter onto 
its surface. In our model, the AXPs experience propeller spindown (with some accretion) 
when Tin > Tco- The characteristic ages Xj of the disk estimated from the propeller torques 
Ni, N3, and iV4 are greater than Tc when > Tco, depending on the value of S*. Figure 3 
shows that the present age estimate is larger than Tc for all of the torques considered here if 
the magnetic dipole field strength on the poles is greater than 6x 10^^ G. If the end of the 



propeller phase is taken to occur when ri. 



as commonly assumed and adopted by 



WCK06, then B^u ^ 0.6 is obtained with ~ 1.3 x 10^ yr = n = = n for Md = 10"^ Mq 
and ~ 1.3 x 10^ yr = n = rg = u for Ma = 10"^ Mq (see Fig. 3). 



In scenarios with these propeller torques A^^i, N^, and A^4, 4U 0142+61 must have spent 
most of its lifetime in the propeller phase with little or no accretion, and it must have en- 
tered its current phase of spin-down with accretion more recently. In the past, all or most 
of the incoming disk matte r would have been thrown o ut of the system as envisaged in the 
original propeller scenario (jlllarionov fc SunyaevI 119751 ). While in the present stage the gas 
disk is supplying the X-ray luminosity of 4U 0142-1-61 by accretion, in the past the X-ray 
luminosity was not supplied not primarily by accretion but by some other mechanism, such 
as dissipation of the magne tic field in the neutron star's surface, as in magnetar models 
(IThompson &: Duncan! Il995l ). or by internal energy dissipation in the neutron star, as ex- 
pected for neutron stars under torque s but without initial cooling or accretion as a dominant 
source of surface thermal luminosity (lAlpar et al.l Il984j : lAlparl I2OOII ) . 



If, on the other hand, there is accretion onto the neutron star over some parts of its 
surface while a fraction of the infalling disk matter esc apes from the stellar magnetosphere 
as outflow or wind (see llUarionov fc Kompaneetslll990l ). the torque is given by 



(9) 



where Mout = M — M^cc by conservation of mass, with the mass accretion rate. Mace = 
LxR*/GM^, expressed in terms of the X-ray luminosity Lx- We have defined the mass 
outflow rate Mout and the accretion rate Mace as positive quantities. When the negative 
propeller torque described by the first term dominates, the net torque A^5 acts to spin down 
the neutron star. Such a spin-down phase allows one to explain the persistent-stage Lx of 
AXPs and SGRs with partial accretion from a fallback disk, as we have done here. The 
first term in equation Q is the braking torque on the neutron star due to mass loss when 
Mout < (jGhoshlll995l ). The second term represents the angular momentum flux onto the 
neutron star due to accreting disk matter. If accretion of matter onto the neutron star has 
already started a switching off the effective propeller mechanism within the recent history 
of 4U 0142+61, we would expect to find commensurate estimates age and magnetic field 



strength based on both N^, and one of the propeller torques A^i, N3, and A^4, whichever 
is an apt description of the actual torque effective in the system's past. In Figure 3, we 
plot the age and polar surface magnetic field strength estimates from A^^i, A^3, A^4, and N^. 
Note that the spin-down torque and the propeller torque A^i yield a common estimate 
{B^ ~ 6 X 10^^ G) for the dipolar component of the magnetic field. Even for the smaller disk 
mass {Md = 10~^ Mq), the torque models N3 and N4 produce very large ages for B^, > 10^^ 
G. The torque represents the most plausible model, because it incorporates accretion with 
spindown and yields reasonable ages, together with the torque Ni, for -Bdipoic ~ 10^^-10^'^ G. 

Using the observed parameters of 4U 0142+61 and torque A^s, the mass inflow rate in 
the disk, for a given magnetic field strength, can be estimated from equation (Q along with 
the outflow rate Mout- We find M ~ 1.9 x 10~" yr-^ and Mout ^ 3 x 10"^^ y^-i 
for 5* = 6 X 10^^ G and Mace ^ 1-6 x 10~^^ M© yr^^ (~ 10^^ g s"^). Note that most of the 
mass lost by the disk accretes onto the neutron star. This estimate implies /2 = 0.85 with 
/i = l. 

The upper limit set on the lifetime of the disk by employing Nr, is much smaller than 
the lifetime values estimated using the propeller torques and A^4. The spin-down age of 
the pulsar estimated with the dipole spin-down formula, which of course can only offer a 
rough comparison with the present model, P*/2P* ~ 10^ yr, is comparable to t^, the inferred 
disk lifetime corresponding to the spin-down torque iVs in the accretion regime (see Fig. 3). 
This age estimate is also consistent with those for the supernova remnants associated with 
the AXPs IE 2259+586 and IE 1841—045. We conclude that an active disk torquing down 
the star for an age of ~ 10^ yr is a tenable hypothesis, most likely with some accretion for 
a significant part of its history up to the present. 

With the assumption that is a good description of the current torque, one can 
estimate the fraction /2 of the mass inflow that is accreted, for a given value of the disk 
inner radius, without having prior knowledge of the mass inflow rate in the disk. Equation 
(9) can be written in the form 

where we have substituted = Using the observed spindown rate and the accretion 

rate inferred from the luminosity (eq. [2]), we find /2 = 0.6 for rin = 10^ cm and /2 = 0.4 for 
rm = = 7 >^ 10^ cm. 



(10) 
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4. DISCUSSION AND CONCLUSIONS 

By combining earlier optical data (HvKKOO, HvKK04; Morii et al. 2005) with recent 
Spitzer infrared data (WCK06), we have shown that the full data set is consistent with the 
luminosity and spectrum expected from a gas disk. The optical radiation is powered to a 
significant extent by viscous dissipation in the inner regions of the disk while the IR radiation 
comes from outer regions where irradiation by the central source is the primary source of 
the disk's local luminosity. 

What is the maximum dipole magnetic moment //* that is compatible with these disk 
model fits to the data? Prom equations (1) and (2), we obtain 

= /r'/V2-^/v,rMV^(GM.)V4. (11) 

Using fi — 0.5 and /2 = 1.0 and employing the maximum rin — 8 x 10^ cm that is compatible 
with our fits, for a minimum Ay — 2.6 we find the maximum possible to be 5.8 x 10^^ 
G cm^ corresponding to a maximum surface dipole field of 1.2 x 10^^ G on the poles of a 
neutron star of mass 1.4 and radius i?* = 10^ cm. A further increase in the upper limit 

for the magnetic field can be obtained with smaller values of /2. The minimum /2 that allows 
a fit within 60% V-band deviation is around 0.3, together with Tin ~ 1.0 x 10^° cm. For these 
extremes of parameter values, we obtain 5^, ^ 3 x 10^"^ G on the pole of the neutron star. 
These estimates should be noted together with the uncertainties due to disk-magnetosphere 
interaction, possible shielding effects, and time-dependent variations as discussed in § 2. 

We note that if the disk inner radius is large, Tin ^ r^o, the neutron star is a fast rotator, 
far from rotational equilibrium with the disk, and likely to be in the strong propeller regime. 
If the simple description of propeller spin-down given by equation (9) is roughly correct, as 
is likely the case, the large moment arm provided by a large Tin, together with mass outflows 
at least comparable to (and probably larger than) the mass accretion rate onto the star will 
yield spindown rates much larger than that observed: 

-MacArf„. (12) 

A disk with Tin ^ 1.0 x 10^° cm, which allows for a magnetic dipole as large as ~ 3 xlO^^ 
G, requires a spindown rate of ~ 10~^° rad s~^. This is much larger than the observed 
~ 1.7 X 10~^^ rad s~^. Here we have a strong indication against values of large 
enough accommodate dipole magnetic fields in the magnetar range. 

As our disk models depend only on the mass infiow rate M inside the disk and the 
irradiation parameter C, it could be argued that the entire mass infiow actually turns into 
an outfiow and none of the mass fiowing through the disk is accreted onto the neutron star. 
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the X-ray luminosity being due to some energy release mechanism other than accretion. The 
observation that the inner disk radius in our best fit is so close to the co-rotation radius argues 
against this. Indeed our fits yield the inner radius required by the optical observations. The 
actual inner radius may be even smaller and closer to r^o, but we do not have the possibility 
of observing such a highly absorbed source in the UV band to check. The rather small value 
of the spin-down rate compared with what one would expect with a substantial mass outflow 
from a larger inner radius also indicates that this is a weak propeller, close to equilibrium. 
It is therefore likely that there is mass accretion. 

Could the disk of our model fits be passive, with no viscous dissipation, no mass inflow, 
and no accretion at all onto the neutron star? Under this hypothesis, the disk luminosity 
would be generated entirely by irradiation by the neutron star's X-ray luminosity, which in 
turn is not due to accretion. It is unlikely that there is no accretion at all, if there is a gas 
disk with the properties indicated by our model fits, which are sensitive to the inner disk 
radius and require temperatures above 10^ K in the inner disk. The disk cannot be neu- 
tral and passive at these temperatures; it has to be ionized, whether it is pure hydrogen or 
contains heavy elem ents. An ionized disk wo uld have viscosity generated by the magnetoro- 



tational instability (IBalbus fc Hawleyi Il99ll ). which will result in mass inflow through the 
disk and interaction between the disk and the neutron star magnetosphere, as well as some 
mass accretion onto the neutron star, depending on the closeness to equilibrium. Ionization 
temperatures for hydrogen disks are of the order of 6000 K. For fallback disks with heavy 
metals (Werner et al. 2007), the ionization temperatures could be as low as ~ 2500 K. In 
hydrogen disk simulations of soft X-ray transients and dwarf novae, two different viscosity 
parameters are employed above and below the ionization temperatures to account for the 
observations of these systems. In these models, the viscosity parameter decreases by a fac- 
tor of 0.01-0.5 as the temperatures decrease to below the ionization temperature, Ti. This 
indicates that even the disk regions with T < T[ remain active and that the critical temper- 
ature below which the entire disk enters a passive state is likely to be much lower than the 
ionization temperatures. Early studies by Gammie (1996) on protoplanetary disks suggested 
the possibility that a weakly ionized disk might become passive in its dense regions in the 
mid-plane while accretion proceeds over its surface. Recent work by Inutsuka & Sano (2005), 
however, indicates the existence of various physical mechanisms that can provide, sustain, 
and homogenize a sufficient degree of ionization for the magnetorotational instability (MRI) 
to work through most regions of a protoplanetary disk even at temperatures as low as 300 
K. According to Inutsuka & Sano, once MHD turbulence starts to prevail in a disk, it seems 
quite difficult to switch it off. The state of turbulent viscosity keeps the temperature and 
other parameters such that the ionization level required for MRI is self-sustained. 



The presence of the disk imposes constraints only on the dipole component of the mag- 
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netic field. The very strong magnetic fields that are employed to explain SGR and AXP 
bursts in the magnetar model may well reside in the higher multipole structure of the neu- 
tron star's surface magnetic field. In this hybrid situation, dissipation of the magnetar 
fields' energy in the neutron star crust may contribute to the surface luminosity along with 
accretion. In this case, Mace will have a value less than that inferred by taking the X-ray 
luminosity to be fully due to accretion, as in equation (2). 

The object 4U 0142+61 is the anomalous X-ray pulsar with the most detailed optical 
and IR observations. It would be of great interest to check and extend the current database 
with further observations, in particular with Spitzer. It is important to realize that magnetar 
fields in higher multipoles and a dipole field with i?^,<10^^ G, to accomodate a disk with 
the equilibrium period range of observed period clustering, can coexist in the neutron star. 
The search for disks in other AXPs has also provided some data in the optical and near - 



infrared that are consistent with gas disk models, as discussed in lErtan fc CaliskanI (120061 ). 
Time-dependent variations in the X-ray luminosity could be smeared out during the X-ray 
reprocessing at the outer disk. On the other hand, variations in the mass inflow rate modify 
the disk emission in the optical bands before the changes have been observed in X-rays, 
because of varying accretion rates. Any observed signatures in the optical bands followed 
by similar variations in X-rays with viscous time delays (minutes to days) would be a clear 
indication of ongoing disk accretion onto the neutron star. At present, 4U 0142+61 seems 
to be the best source with which to test this idea. 

Finally, we comment on the passive dust-disk model proposed by Wang et al. (2006). 
Their detection of 4U 0142+61 in the mid-infrared has firmly indicated, for the first time, 
the presence of a disk around an isolated neutron star. They fitted the data with a two- 
component model, using a power law to fit the optical and near IR data, which they ascribe 
to magnetospheric emission. This leaves only their recent discovery, the Spitzer detections in 
the mid-IR band, to be explained by a disk. Hence, while the evidence in that band clearly 
points to a disk, the association of only this narrow-band emission with the disk constrains it 
to be a dust disk beyond the magnetosphere. Questions arising from this model concern why 
this disk did not generate some viscosity to extend into an active gaseous disk, how it cooled 
and remained cold and confined, and whether it is stable against the radiation pressure of 
the pulsar. With an inner disk radius r^^ ~ 2ric, where rj^ is the light cylinder's radius, 
such a passive dust disk could be stable against the radiation pressure of the magnetospheric 
pulsar emission (Ek§i & Alpar 2005). This stability is valid for ric ^ r^^ < 2ric for orthogonal 
rotators and for even wider range for non-orthogonal rotators according to the investigation 
of Ek§i & Alpar (2005), the stability issue should be kept in mind for dust disks beyond the 
light cylinder. 
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Table 1: Model Temperatures and Radii Corresponding to Different Optical and Infrared 
Bands. 



Band 


Tbb (K) 


R (cm) 




B 


6516 


7.0 X 10^ 


0.45 


V 


5263 


1.0 X 10^° 


0.3 


R 


U54: 


1.4 X 10^° 


0.30 


I 


3585 


2.2 X 10^° 


0.14 


J 


2377 


4.8 X 10^° 


0.07 


H 


1779 


7.9 X 10^° 


0.04 


Ks 


1324 


1.4 X IQi^ 


0.02 


4.5 fim 


644 


5.9 X IQii 


0.005 


8 iim 


362 


1.9 X 10^2 


0.002 



Note. — 



The rightmost column shows the ratio of viscous dissipation rate to irradiaton flux 
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Fig. 1. — Energy flux data and irradiated disk model values for 4U 0142+61 in the optical 
and infrared bands (B, V, R, I, J, H, K^, 4.5 //m and 8 //m). Circles are the observed 
(absorbed) data (taken from HvKKOO [V, R, I], HvKK04 [B, K,], Morii et al. 2005 [H, J], 
and WCK 2006 [4.5 /xm and 8 /^m]), and triangles are data dereddened with Ay — 3.5. 
Squares are the irradiated-disk model energy flux values (see § 2). 
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Fig. 2. — Same as Fig. 1, but with dereddennig of Ay — 2.6 (triangles). 
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Fig. 3. — Disk lifetime as a function of stellar magnetic field strength for different propeller 
torques. Squares show the disk lifetime values estimated from the spin-down torque in 
the accretion regime. The left vertical axis shows the estimated lifetime range for a disk of 
mass Md = 10~^ M©. The vertical axis at right is for Ma = 10~^ Mq. 



